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ORBIT IHU ALINEMENT IKTERPRETATION 
OF ONBOARD DISPLAY DATA 

By Roger Corson, HcDonnell Douglas 
Technical Services, Co., Inc. 


1.0 ^UHUAjaX 

This document Investigates the Space Shuttle Inertial measurement unit (IMU; 
allnement algorithm, to determine the most Important allnement starpalr selection 
criterion. Three crew-displayed parameters are considered: (1) the results of 

the separation angle difference (SAD) check for each starpalr, (2) the separa- 
tion angle of each starpalr, and (3) the age of each star measurement. It was 
determined that the SAD for each pair cannot be used to predict the IMU 
allnement accuracy. If the age of each star measurement Is less than approxi- 
mately 30 minutes, time Is a relatively unimportant factor and the most Impor- 
tant allnement pair selection criterion Is the starpalr separation angle. 
Therefore, wt?en there are three available allnement starpalrs and all measure- 
ments were taken within the last jO minutes, the pair with the separation angle 
closest to 90 degrees should be selected for IMU allnement. 
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2.0 INTRODUCTION 


The purpose of this paper Is to deteralne how the Shuttle IMU allnemerit accuracy 
can be predicted, using the star neasurement data that are available to the crew 
via the star tracker/crew optical allnement sight (COAS) control monitor cath- 
ode*ray tube (CRT) display. The allnement algorltnm will be analyzed In order 
to form conclusions relating to the crew and/or ground personnel Interpretation 
of displayed measurement data for the purpose of allnement starpalr selection. 
Several allnement-related characteristics of the current version of the star 
tracker software will be discussed briefly. 
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3.0 piscgssiQN 

Alinement of the Space Shuttle IMU platfoms to a dealred reference orientation 
Is a critical navigation procedure during onorblt operations. Nominally, each 
Inertial platform Is alined to Its desired orientation by first determining Its 
present orientation and then repositioning It to the desired orientation. The 
present orientation of each IMU platform Is determined by taking directional 
measurements on a pair of stars. The star measurements are expressed as 
unit llne*of-slght (LOS) vectors In a coordinate system that Is fixed relative 
to the IMU platforms. Nominally, these star measurements will be taken using 
the Shuttle star trackers. These electro-optical tracking Instruments are 
used because of their accuracy and their ability to be used In an automated 
system. 

During the first two orbital flight test (OFT) missions, star tracker measure- 
ments will be taken manually by the crew. The majority of this discussion, how- 
ever, concerns the automatic op'>ratlon of the star trackers that will be 
employed after Initial OFT missions. Under computer software control, each of 
the two trackers will acquire and track stars of opportunity. Measurement data 
for each star tracked are evaluated for validity by the software and then saved 
as a LOS vector In a star-sighting table. This sighting table Is sized to store 
data on a maximum of three stars. After three stars have been saved, with each 
additional star tracked, the software determines which of the four stars Is of 
least value for alinement computations and discards that star, leaving Just 
three In the table. 

Data for only two stars are needed for computing the present orientation of each 
IMU platform. Because there are Inherent errors In the measurements, there 
exists an optimum starpalr geometry which will result In the minimum IMU 
alinement error. This alinement error relation to the pair geometry Is 
discussed In sections 3.1 and ?,.2. Since stars are acquired In a more or less 
random fashion (rather than according to seme optimum, deterministic pair selec- 
tion criteria as Is done In the ma..ual sighting mode), the stars appearing In the 
sighting table will have random and unpredictable pair geometries. By 
Increasing the number of stars saved from two to three, the number of available 
alinement pairs :ncreases from one to three, thereby Increasing the probability 
(at any one time) of having a pair with an acceptable geometry. 

For each star that Is saved In the sighting table, several related parameters 
are displayed to the crew on the star tracker/COAS control monitor CRT display. 
The identity number of each star that has been tracked and saved In the sighting 
table is displayed. Other Information displayed Includes the time since each 
star was sighted, and the separation angle of each of the three pairs. One 
final display Item Is an error parameter for each of the three starpalrs. The 
SAD error (ERR on the display) 1s the absolute difference between the separation 
of the two measureo LOS vectors and the actual or known separation between the 
two stars. When the software determines that there Is at least one pair of 
stars existing In the sighting table that has an acceptable geometry for 
alinement, the software will select the two stars having the beat pair geometry 
and flag these two stars on the crew display as Its chosen pair for alinement 
computations. The crew may then execute the IMU alinement using the software 
selected stars. Alternatively, a crew or ground personnel review of the star 
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tracKer display parameters may result In the selection of a different pair of 
stars to ensure an optimum IHU allnement. The purpose of this discussion Is te 
analyze the allnement algorithm in order to form conclusions relating to the In- 
terpretation of erew displayed star data for use In IHU allnements. 


3.1 IHU ALINEHENT ERROR 

The IHU allnement error relations will be derived geometrically by Inspection. 
Each IHU will have some random orientation with respect to the measured 
starpalr. The allnement errors, therefore, will be arbitrarily referenced rela- 
tive to the coordinate system shown la flgu.’e 1. The X-Y plane Is defined by 
the actual starpalr LOS vectors S and T and will henceforth be referred to 
as the pal” plane. The X-axis bisects the angle between the starpalr LOS 
vectors. The Z-axls la perpendicular to the pair plane and the Y-axls com- 
pletes the right angle triad. This orientation was selected to simplify the 
derivation of the error equations. 

The measured star LOS vectors S and f are perturbed from the actual by small 
errors. The errors In the measured star LOS vectors are resolved Into compo- 
nents parallel (II ) and perpendicular (1) to the pair plane. When an allnement 
Is performed, the computed or measured orientation will be In error as a result 
of the sighting errors. The Shuttle allnement method assumes that the measured 
starpalr lies In the same plane as the antual starpalr (ref. 1). This method 
further assumes that the measured and actual directions of one of the pair are 
coincident. This direction vector will be referred to as the principal LOS of 
the pair. In this analysis the vector S will be the principal LOS. The 
allnement error, therefore, for this example Is measured between the plane 
defined by S and T and having a reference direction S and the plane 
defined by S and T with the corresponding in-plane reference S. An equa- 
tion for each component of the allnement error vector ^ will be derived. 

The error about the X-axls, Is a result of the difference between Sj^ and 
T|, the out-of-plane errors. If both the In-plane and out-of-plane errors 
are small, then the tamgent of Is equivalent to the dlf.^erence between 

Sj^ and Tj^ divided by the distance between the end points of S and T. 


4'jf w tan Vj 


Sj - 

2 sin ^/2 


6 


5 starpalr separation angle 
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The error about the Y-axis, iy, Is a result of the average nagnltude of Sj^ and 
T|. Assualng that the errors^are snail, the tangent of ♦ , therefore. Is 
equivalent to the average of 8| and divided by the distance between the 
end point of either S or T and the Y-axls. 


♦ y « tan ♦y 


1/2(S| ♦ Tj^) 
ooa */2 


The error about the Z-axls, ♦j. Is slaply a result of the component of the error 
In S that lies In the pair plane. 


♦z - S|| 


In mumnary, the components of the allnenent error vector are 
♦x » - Tx)/sln */2 

<ty » 1/2(S| Ti>/coa */2 

♦z * S|| 

Assume all sighting errors are independent Qaunslan with zero 
standard deviations. 


(la) 

(lb) 

(lc) 

means and equal 


Oq = 10 error In S||, Sj^, T||, Tj^ 


The one- 3 igma error In (Sj^ -Tj^ ) or (S| + Tj^) Is equal to The errors In 

the mlaalinoments <t>^, <t>y, oil have zero moan and standard oevlatlon 

Ojj Oy, Oj respectively. 


1/2 

Oq v^i/sln ^/2 

(2a) 

1/2 

0 ^^ 2/003 '5/2 

(2b) 

°o 


(2ft) 
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A tot;_I allneaent arror Indicator will ba dafinad as a r':ot sua squara (RSS) of 
aaoh of tha ooaponants and 


♦ « ♦ ♦y*^ ♦ 

Assuming this dafinition, it can ba shown that tha root mean squara (RMS) of ^ 
is equivalent to an RSS of tha standard deviations of each of tha components 
Ojj, Oy, and 

RMS<> 5 <*) 2 \ Oy2 ^ 0^2 

'»» « Oq Jy * 2 csc^6 (3) 


It can also be shown that the standard deviation of tha total sighting error ('^) 
is related to each axis standard deviation (0^) by: 


o . -J^v'2 


The total RMS alinement error expressed as a function of the total system 
sighting error and the stnrpair separation, therefore, is: 


u) s o 



♦ csc^6 


(4) 


The RMS indicator ia equivalent to an RSS of the mean and standard deviation of 
the total alinement error. This equation for the RMS alinement error (eq. (4)) 
agrees with that derived by R. P. O'Donnell of the Massachusetts Institute of 
Technology (MIT) in reference 2. In figure 2 the total RMS alinement error, o), 
is plotted as a function of both ^ and In figure 3 the standard deviation 

of each component of the alinement error la plotted separately. 


The alinement error equation (eq. 4) was derived, assuming that both star 
sightings were taken at the same time and, thus, does not Include the effects of 
errors duo to IMU drift. If IMU drift error Is added to each star measurement, 
the alinement error equation becomes quite complex. Ro'irlting equation (1) to 
include drift errors yields 

(J)^ = 1/2 (Sx -if )/sin ^/2 
<iiy = 1/2 (Si Ti)/cos ^/2 
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Alinenent starpair separation angle, 6 (deg) 
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wher« the barred errors are defined: 

Si 2 ^i ♦ drift error 
Ti - Ti ♦ drift error 
S|| 2 S|| ♦ arlft error 

The drift errors are a function of the following parameters: 

tg 2 time since the sighting of S 

t^ S time since the sighting of T 

djj 2 IMU drift rate about X-axls 

dy = IMU drift rate about Y-axls 

d_ 2 IMU drift rate about Z-axls 

There fore , 

Si . S| ♦ tg d^ sin ^/2 ♦ tg dy cos */2 

fi s Ti - t^ d^ sin ®/2 ♦ t^ dy cos ^/2 

^11 ' ^11 ^ 

These equations are derived from the geometry as illustrated in figure 1 . The 
three perpendicular components of the alinement error, therefore, are: 

« 

(J)^ = 1/2 Si - Ti ♦ d^ sin ^/2(tg ♦ t^) ♦ dy cos */2(tg - t^.) /sin */2 

^y = 1/2 Si ♦ Ti ■*. d^ sin */2(tg - t^) dy cos */2(tg ♦ t^.) /cos ^/2 

■ ^11 ^3 ^Z 
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Now, ooaputinK the standard deviacion of each coaponent yields (again, the nean:i 
of ^y, and ^2 shown to be r.ero): 

Ojf M 1/2 y/ 20 ♦ ( tg^ ♦ t^^ - 2tgt^ oos 

Oy w 1/2 ^20^2 ♦ (tg2 ♦ tt^ ♦ 2tgtt eoe 

Oj . ^ ^ ♦ t o/ 

* V O B d 

Recall again that the standard deviation of the total sighting error Is related 
to each axis standard deviation by: 

o . a^y/2 


/sin ^/2 
/cos ^/2 


The total RMS allnement error expressed as a function of the total system 
sighting error, the pair separatloi, and the age of each sighting, therefore, 
is: 


'!) 


L 


12 ♦ cst"^6 1/2 


r 1 

— [tg2 ♦ (tg^ ♦ tt^)CSC^6 - 2tgtt COt ^6 I 


(5) 


3.2 STAR TABLE CRT DISPLAY DATA 

If both star sightings were taken fairly recently, the allnement error contribu- 
tion resulting from IMU drift will be small. Tnls Is because the ratio 
(eq. (5)) Is very small for the Shuttle system (tt5 x 10"^/sec). Furthermore, 
for the majority of the time the age of one of the sightings In each pair will 
be zero and the age of the other star will be less than 30 minutes. Simulations 
cf the automatic acquisition mode have demonstrated that sighting data older 
than 30 minutes rarely appear in the table and most of the star data is less 
than 15 minutes old. Consider a pair of alinement stars with a separation angle 
of 90 degrees - the age of the most recent sighting is zero, and the other star 
neasurement Is either zero, 15, or 30 minutes old. The 30-minute case yields an 
'M.'! allnement error only 0.27a greater than the optimum or the zero-minute case 
(eq. (5)). This is approximately 20 arc seconds for the Shuttle system. In the 
15-minute case, the allnement error as predicted by equation (5) is only O.OTo 
or approximately 5 arc seconds greater than the optimum case. For fairly recent 
sighting data, therefore, the pair separation an.;le affects allnement accuracy 
to a greater degree than sighting age does. (fig. 2). 

Before concluding this discussion of the effect of IMU drift on alinement error, 
one final important point should be noted. The time since the measurement of 
the principal LOS of the starpalr, t^,, appears in one term In equation (5) all 
by Itself. The alinement error, therefore, is more strongly related to the IMU 
drift error in the principal star LOS. Since the principal LOS is favored by 
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the allnement algorithm, this LOS should be the one that was most recently 
sighted. This procedure can result In an Improvement In allnement accuracy of 
up to 20 percent . 

The most obvious property of the curve In figure 2 demonstrates that the 
allnement error Is significantly degraded for star separation angles close to 
zero or 180 degrees. It Is desirable, therefore, to select allnement starpalrs 
with near-right angle sepa.'atlons In order to guarantee minimum allnement error. 
Crew selection of an allnement pair should be driven by the separation angle 
displayed for each pair. The current version of the star tracker subsystem 
operating program (SOP) (ref. 3), however, already performs this selection based 
on pair separation, thereby relieving the irew of any decisive action if the 
data have been taken recently. The crew should, however, verify that all data 
In the sighting table are less than 30 minutes old before accepting the software 
selected pair. 

There Is one final Item that Is a candidate parameter for use In predicting 
allnement accuracy. This parameter, the SAD, Is the absolute difference between 
the measured pair separation and the actual or known pair separation. The SAD 
error that is displayed to the crew Is defined: 


ER t s 



(4) 


The SAD error Is clearly a function of only the i.i-plane errors. If the SAD 
gave any Indication of the magnitudes of the lr.-|*lane errors, it would still bj 
relatively invaluable because for angles not cloae to 90 degrees, the out-of- 
plane errors are the greater contributors to the allnement error (fig. 3). How- 
ever, assuming tnat both stars of the pair were correctly Identified, closer In- 
spection of the equation for ERR reveals that knowledge of the SAD gives no clue 
to the magnitude of the In-plane allnement error. Solving equation (4) for 
the In-plane allnement error 


'^z ■ ^11 = “^11 t 


Since T|| Is a random variable, a better estimate of the In-plane allnement 
erre*, cannot be determined using the known or measured value of ERR. 


The displayed SAD error cannot tc used to predict allnement error when both com- 
ponents of the starpalr have been correctly Identll’icd. Th», SAD check Is 
performed in the onboard Shuttle software to detect failure of the star tracker 
to acquire the desired (or so ft ware- selected star). The use of measurement data 
on the wrong star would result in very large allnement errors. To protect 
against this possibility, the star tracker software checks the SAD of each new 
star tracked and the previously tracked star (If available) for off-nominal 
values of ERR. If the value of ERR is very much greater than Its nominal value 
^*^ERR " least one of the stars in the pair was fal"ely Identified. 

If the test fails, the new star data is not saved in the star sighting table. 
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Putheraore, if the test Pells and the previously sighted star Is the only one in 
the sighting table, then It Is deleted also. If the SAD check Is passed, the 
new data Is stored In the table. The software assumes that If a pair of stars 
passes this check, there la a high degree of certainty that both stars ha\ » been 
correctly identified. For this reason, each star Is paired with onl/ one other 
star for tl»e purpose of performing the SAD check. It Is erroneous to assume 
that the software SAD check provides absolute certainty of the star Identities. 
This Is because the SAD check Is only sensitive to random errors that lie In the 
star plane (ref. U). As stated previously. If the table Is al'*eadv full, the 
software determines which star In the table is of least value for .tlln<riiient 
computations, and deletes that star to make room for the new star. Any one of 
the three stars In the table can be replaced by the new star. The SAD for each 
pair of stars In the sighting table Is displayed on the CRT. When a new star Is 
placed Into a full table, the set of SAD angles may be recon .'Igu red In such a 
way that star pairs previously tested by the software SAD check no longer exist. 
By this software sighting table maintenance scheme, therefore. It Is possible 
that one, two, or all three of the displayed SAD angles have not been tested by 
the software against the allowable tolerance. Figure 4 Illustrates how multiple 
pairs with common LOS vectors produce Intersecting tolerance limits about each 
star LOS. Out-of-plane errors In a particular pair of stars can be seen by the 
SAD errors of other pairs with a single common LOS vector. It would be advanta- 
geous, therefore, for either the crew or ground personnel to check the displayed 
SAD angles for all pairs In the table that have not been checked by the soft- 
ware. This procedure also may Isolate out-of-plane errors that were not 
detected by the software SAD check. If the SAD check Is not done for the pur- 
pose of detecting out-of-plane errors, the tolerance windows will be bounded In 
the direction perpendicular to the pair plane by only the estimated/measured po- 
sition check tolerance. This tolerance. Illustrated by the circles In figure 4, 
however. Is very coarse compared to the SAD error tolerance. The estimated/ 
mea.iured position check tolerance Is currently set at 0.5 degrees for offset 
acqui:?ltlon and 1.8 degrees for automatic acquisition. 


SAD limits 



Figure 4.- Three pair SAD tolerance windows. 
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4 . 0 CONCLUS IONS AN D _R|COWENDAT WNS 

During th« first two OFT missions, IMU allnements will b« done manually. The 
nominal mode of operation will be to maneuver to an attitude In order to 
simultaneously track one star In each star trackw. The starpalr separation 
angle will always be very close to 90 degrees and the age of both sightings will 
always be near zero. Therefore, the Information presented here la not reltfvant 
to the Initial OFT missions. During later missions however, when IMU allnnments 
are done using the automatic acquisition mode, the following guideline Is 
recommended whenever there la data for three stars stored In the sighting table: 
If all sighting data were taken within the last 30 minutes, the starpalr separa- 
tion angle Is the moat Important allnement pair selection criterion. Pair se- 
lection based on separation Is already implemented In software, and the crew 
should not override the software decision by selecting a different pair. If any 
single star In the table Is older than 30 minutes, optimum pair selection be- 
comes complsx because of the Interrelation of allnement accuracy, time, and pair 
separation. Development of a simple allnement accuracy Indicator based on time 
and pair separation angle (for use In pair seletlon) Is a subject of future 
study. 

Two allnement-related software characteristics that were discussed are mentioned 
again In summary. First of all. If the most recently sighted star Ic used as 
the principal LOS, this can result In a considerable improvement in the 
allnement accuracy over the converse case. This Implementation Is reflectec by 
the IBM Corp. In their current Interpretation of the star tracker SOP. Sec- 
/•Mdly, when three stars have been saved In the sighting table, a greater degree 
of certainty of the Identity of each star sighted can be galncu by having the 
crew or ground personnel monitor the SAD error for all .'•airs that have not been 
checked by the software. If the displayed value of ERR for any pair la greater 
than the I-load SAD tolerance (TOL10), then at least one star In the pair has 
been misldentif led . 
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